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Abstract 

There remains a useful relation between diffusion and mobility for a 
Langevin particle in a periodic medium subject to nonconservative forces. 
The usual fluctuation-dissipation relation easily gets modified and the 
mobility matrix is no longer proportional to the diffusion matrix, with 
a correction term depending explicitly on the (nonequilibrium) forces. 
We discuss this correction by considering various simple examples and 
we visualize the various dependencies on the applied forcing and on the 
time by means of simulations. For example, in all cases the diffusion 
depends on the external forcing more strongly than does the mobility. We 
also give an explicit decomposition of the symmetrized mobility matrix as 
the difference between two positive matrices, one involving the diffusion 
matrix, the other force-force correlations. 

1 Introduction 

The relation between internal fluctuations in a system on the one hand and 
its susceptibility under an applied external field on the other hand is of great 
practical and theoretical interest. The most immediate example concerns the 
relation between diffusion and mobility and remains important in the develop- 
ment of nonequilibrium statistical mechanics. The first studies go back to the 
works of [25l [26] and of [7], (see also [19]). They are early examples of the 
fluctuation-dissipation theorem, itself being a cornerstone of linear response 
theory for equilibrium systems, see [H] and 

When moving to nonequilibrium systems, the fluctuation-dissipation relation 
(FDR) is typically violated and there is no a priori reason why the system's 
mobility can be simply obtained from its diffusivity. All the same, mobility 
remains crucially important for the discussion of transport properties under 
nonequilibrium conditions |101 124j . That is why we want to get explicit infor- 
mation about the corrections to the FDR, a subject that has been considered 



from many different sides. We mention few. One approach has been to as- 
sociate an effective temperature to the nonequihbrium system, for example in 
ageing regimes, as e.g. in [l] and [21]. Another approach is to modify the FDR 
inserting the nonequihbrium stationary density, or going to a "moving" frame 
of reference using the probabihty current or local velocity [S] , somewhat in the 
spirit of the approach of [5] , see also [TTJ [TH [13 US] . Specific models of inter- 
acting particles have for example been treated, also in mathematically rigorous 
ways, by [H [mini HI]. More recent experimental work includes [TUl [201 [5T] . 
Here we concentrate on the nonequihbrium regime where the Sutherland-Einstein 
relation gets modified because of nonconservative forces, starting from the iner- 
tial regime. As will become clear, not much effort is needed for finding examples 
where these corrections are visible and detectable. Recent work by [T31 [H] is 
similar to ours but in a different setting. The present paper gives the correc- 
tions following recent work by [3j |4j , which enables us to zoom in on particular 
dependencies, such as on environment and driving parameters. 
The main results of the present analysis are the explicit relations ([7| and ( 10 1 be- 
tween the mobility and the diffusion constant, and their numerical exploration. 
In particular, we typically find that for the models treated here the diffusion 
constant depends on the external (nonequihbrium) forcing in a stronger way 
than the mobility. The mobility is bounded from above by the diffusion, as we 
derive in an exact bound ( 20 ) . When the forcing only depends on the position, 
then the correction to the Sutherland-Einstein relation is second order. 



The plan of the paper is as follows. Section [2] contains the general relation 
between mobility and diffusion matrix for Langevin dynamics in a magnetic field 
plus other general external forces. The correction to the Sutherland-Einstein 
relation is in terms of a correlation function between the particle velocity and 
the applied forces. As simplest and standard illustration Section [3] reminds us of 
the theoretical framework of the Sutherland-Einstein relation. It is the famous 
proportionality of the mobility (in linear regime) with the diffusion constant 
(of the original system). We see from these textbook examples that our ap- 
proach is suitable. Then restarts the nonequihbrium analysis. Section |4] treats 
nonequihbrium fluids with uniform temperature. The nonequihbrium condition 
is imposed by stirring the fluid (rotational or nonconservative forces) . Here our 
motivation is primarily to illustrate the existing exact formulae and to explore 
more the role and the influence of the various parameters. All involved observ- 
ables can be measured under the nonequihbrium averaging, and vice versa, the 
formulae in principle enable to learn about unknown driving forces and parame- 
ters from the correction to the Sutherland-Einstein relation. Section [5] explores 
the symmetrized mobility; the diagonal elements allow some general bound. Fi- 
nally, the Appendix contains some more technical points dealing with analysis 
and numerics. 
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2 A general mobility-diffusion relation 



Consider a particle of mass m, which diffuses in a heat bath (e.g. some fluid) 
in according to the Langevin dynamics 

n = = V, (f) 

mvt = w ^*^"!^^ = F(ft, ut) - 7rn7;i + y^2mjT^t 

Here is the position of the particle at time t, and Vt is its velocity. The 
particle is passive and undergoes the influence of a heat bath in thermal equi- 
librium at temperature T. In the usual weak coupling regime the particle suf- 
fers friction with coefficient 7 and random collisions, here represented by the 
vector of standard Gaussian white noises. This means that each of the com- 
ponents has a Gaussian distribution, with mean zero (^t,i) = and covariance 
{^t,iCs,j) = Si,jS{t—s). Throughout this text we set Boltzmann's constant to one. 
Furthermore there is an external force F working on the particle. Throughout 
this text, we restrict ourselves to forces that depend periodically on the posi- 
tion r. We make this restriction to have a form of translation invariance in the 
system. In other words, we consider time-scales over which there is no confining 
potential. The force F can depend on the velocity as for a magnetic field when 
the particle is charged. 

When given an initial (probability) density /i(r, v) of independent such par- 
ticles at time zero, it changes in time according to the Fokker-Planck equation 



— - = ~v ■ VrfJ-t - Vi 
ot 



F — \ 77" ^ 

Mt V^^t 



(2) 



m J m 

There are two important quantities that characterize the transport behavior of 
such a system. The first quantity is the diffusion (matrix) function D{t), which 
is defined as ^ 

Dioit) = ^((^t - {ft - ro)j^ 

The subscripts denote the components of the corresponding vectors, and the 
right-hand side is a truncated correlation function: for observables A and B 

A;b) = (^Ab)-(^a)(^b) (3) 

For diffusive systems, i.e., for the systems described above, this diffusion func- 
tion is expected to have a large time limit, called the diffusion matrix 

Dij = lim Dij{t) 

f 00 

In words, the (co)variance of the displacement of the particle is linear with time 
(for large times t 3> I/7) with slope given by the diffusion constant. Some anal- 
ysis is found in Appendix [Aj Secondly, there is the mobility (matrix) function 
M{t), defined as follows: we add to the dynamics in ([T]) a constant (but small) 
force /, replacing F{f, v) — > F{f, v) + f. The mobility then measures the change 
in the expected displacement of the particle: 



M,-(0 = J|.((n-rU) 
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where the superscript / wants the average to be taken in the dynamics with the 
extra force /. Again, this function is supposed to have a large-time hmit, which 
is called the mobility: 

Mij = lim M„(i) 

f oo 

i.e., the mobility is the linear change in the stationary velocity by the addition 
of a small constant force. 

In the special case of detailed balance dynamics, these two quantities are 
related by the Sutherland-Einstein relation, 



M^j - ^A, (4) 



which is an instance of the more general fluctuation-dissipation theorem. 
When the system is not in equilibrium, mobility and diffusion constants are no 
longer proportional. Additional terms show up: for the dynamics defined in ([T]) 
we get 

M^At) - ^A, W + - roh; (5) 

where the vector ^> is explicitly given by 

* = m{vt - vq) - / dsF{fs,Vs) (6) 
Jo 

The relation between the mobility and diffusion functions is thus modified with 
respect to the large time limits in equilibrium systems by the addition of an 
extra term. This term is the (truncated) correlation function between the dis- 
placement and the functional This functional is explicitly expressed in terms 
of the velocity of the particle of interest and the forces that act on it. In the 
limit of large times, the relation ([5| simplifies somewhat: in this limit we get 

Mi,- = —Dii + lim — - — ((ft - fQ)t;^i) (7) 

where the vector $ is explicitly given by 

$ = - / dsFifs,Vs) (8) 
Jo 

Indeed, the term with the correlation between the displacement and the change 
of velocity does not contribute in the large-time limit: 

lim I (in - fo)i; {vt - vo)j) = (9) 
For an argument see Appendix [Aj In Section [s] we further rewrite ([t]) to obtain 

where each term represents a symmetric matrix. Under conservative forces the 
second and the third term sum to give the first term on the right-hand side. 
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Otherwise, in diffusive nonequilibrium tfie correction term to tlie Sutherland- 
Einstein relation is nonzero with the symmetrized mobility matrix as the explicit 
difference between a diffusion-related matrix and the force-force covariance ma- 
trix, as claimed in the last line of the abstract. A further bound on the sym- 
metrized mobility is added in Section [5j 

Our modified Sutherland-Einstein relation has the advantage of being ex- 
plicit in the dynamical variables. Still the formula is sufficiently complicated 
and some aspects of it require careful examination. We can rewrite formula ([t]) 
into 



1 • 1 /■* 

Mji — —Dji — lim / ds 

^ T ^ i^oo 27mT io 



{^^^■,F,irs,^s)) (11) 



As we see, the correction to the equilibrium mobility-diffusion relation is mea- 
sured by a space-time correlation between applied forcing and displacement. 
One simplification is to look close-to-equilibrium. There we see by time-reversal 
symmetry applied to the reference equilibrium that the correction is only quadratic 
in the applied forcing when F only depends on r. In other words, the devia- 
tions with respect to the Sutherland-Einstein relation are then second order. 
Going further from equilibrium also the inverse question becomes interesting, 
to characterize the nonequiilibrium forcing that produces a given experimen- 
tally determined mobility-diffusion relation. We hope to see in future work that 



formulae like ( 11 ) are also useful to learn about the driving conditions from mea- 



surements of both mobility and diffusion. Numerically, of course the right-hand 



side of (11 1 is easier to determine than the left-hand side. 



Usually, as also in what follows, one considers the diffusion and the mobility 
of a single particle. It is however relevant and also possible to include interac- 
tions with other particles. The equations ([I]) must be changed for the force to 
include this dependence on the state of the other particles, but the main results 
remain unaffected. After all, one can simply redo all calculations in larger di- 
mensions (3A^ for N particles) and consider independent white noises acting on 
all 3N-components. Furthermore, the force F can be time-dependent and for- 
mulae ([?])- ([To I remain the same but where F has further time-dependencies. 



be it explicit or be it because of the time-dependent state of other particles. 



The explicit relations above have been derived in [4j following a path-integral 
formalism. In the present paper we provide a number of examples to discuss 
this relation between mobility and diffusion out of equilibrium. For each of 
these examples, simulations have been done to help visualization, and to inform 
us about the relative importance of the various parameters and terms in ([t]). 
On the other hand, the simulations raise new questions, not all of which we 
are able to answer. The Appendix adds more theoretical considerations to 
the relation between diffusion and mobility, and gives information about the 
simulation method. 
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3 Classical equilibrium theory 



Pure diffusion 

The simplest form of the Langevin equation is the case where all the external 
forces are set to zero 



mvt 



(12) 



The mobility and diffusion can be computed explicitly, by first integrating the 
Langevin equation (adding a force /), 



vt = Woe ' H 

7m 



Jo 



and then using the properties of the Gaussian white noise. Because ( ) =0, 
we immediately get the mobility by again integrating over time, 



1 - e-'^* 



1 



7TO j^mt 



As {£,i,t^j,s) — Si,jd{t — s) the velocity-velocity correlations equal (for / = 0) 



e-7(t+s) ^_ ^MZ^g-7|t-s| 



Simply integrating over time then gives the diffusion, 



A,(0 



1 



TS,. 



1 - e-''* 



1 



1 - e-''* 



7m j^mt 



Taking the limit of large times, we see that 



7m i-i-oo T 



T 



A periodic potential 

When a force is added that derives from a periodic potential U, 

n = Vt 

mvt = -^rU{rt) - jmvt + y/^rnrff 
then, the general expression ^ gives 

M^J{t) = lA,(<) + ^(rt-ro).;(«t-i/o)j) 



(13) 



(14) 



The second term on the right-hand side tends to zero for large times, as said 
before, as does the third term: 



lim 

t— f 00 



\ j ds((n-rQ\-^r,U{r)) (15) 
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yielding the standard Mij — Dij/T; see also Appendix [A] Of course the mobil- 
ity and diffusion are no longer equal to 6ij/{"fm). We get a summary in Fig. [T] 
Observe for example that the mobility decreases with the amplitude of the con- 
servative force, which is logical, since the particle needs to escape potential wells 
to have a non-zero velocity. 



In fact, in [S] an explicit expression was derived for the diffusion for over- 
damped Langevin equations with periodic potentials. The overdamped Langevin 
equation is what one gets when the friction is high and the mass of the particle 
is small, so that inertial effects can be ignored. Formally, one lets 7 — ?> oo and 
TO — >■ 0, while keeping x = constant and finite. The result is that one can 



put mvt to zero in (13). The resulting Langevin equation can then be written 



solely in terms of the position: 

The formula of ,8] for one-dimensional diffusion is given by 



D = 



(16) 



where R is the period of the potential, and the integrals in the denominator are 
over one period. We give a new and short proof of ( 16 ) in AppendixjB] In Fig. [l] 
we have plotted this explicit expression, for a potential U = Acos{x) for various 
values of A. This curve corresponds nicely with the mobility and diffusion we 
got from simulations, where we took 7 = 10 and m = 0.1. 
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Figure 1: Diffusion in potential U{x) — A cosx. In {a) A = 1 and the evolution 
through time is shown. Here T — j — m — 1. C{t) here represents the last 
two terms on the right-hand side of (14 1. In (b) the long time limits of the 
mobility and diffusion are shown for values of A ranging from to 2. Ml and 
Dl are the mobility and diffusion for the case T ~ j ^ m ~ 1, while M2 and 
D2 correspond to T = 1, 7 = 10, to = 0.1. Finally D3 is the diffusion that is 
explicitly calculated by (16). (Online version in color.) 
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A magnetic field 

A first choice of a velocity dependent force in ([T]) could be a magnetic field. 
If we add to ( [l2| a magnetic field, say a homogeneous magnetic field in the 
z-direction, — {xt, yt, zt), B — (0, 0, B), with unit electric charge, 



mvx,t = Bvy,t - jmvx^t + y^2m^T^x:,t 
miiy^t = -Bva:,t - imvy^t + V^rnrff $,y^t 



(17) 



then, ([t]) immediately gives 
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Further explicit calculations show that 
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The diagonal elements of M and D are proportional, but not so for the off- 
diagonal elements. Therefore, adding a magnetic field to an otherwise equilib- 
rium dynamics breaks the equilibrium Einstein relation. This may be counterin- 
tuitive, because adding a magnetic field leaves the equilibrium Boltzmann distri- 
bution intact. In Fig.[2]the results of simulations are shown. We took the initial 
position equal to zero, and the velocity is Maxwellian p{v) oc exp(— mw^/2T), 
which is the stationary velocity distribution. The K(t) show the right-hand side 
of (5). As one can see, they coincide nicely with the mobility (left-hand side 
of (5|). Note that the oscillations of the mobility and diffusion have a period 
2Tim/{BT). 
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Figure 2: Diffusion with a magnetic field in the z-direction. Here T — j — m — 
1, while B = 5 for Figure (a) and B ^ 20 for Figure (b). (Online version in 
color.) 
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4 Nonconservative forces 



Recall that we work in unbounded spaces with periodic forces. In that sense we 
consider diffusion on the torus. A conservative force is then a force for which 
every contour integral that starts and ends in equivalent points equals zero. 
Equivalent points have the same coordinates modulo the period of the force. 
If not, we call the force nonconservative. Practically, we distinguish two kinds 
of nonconservative forces. First of all there are forces that cannot be derived 
from a potential (or equivalently, that have a nonzero curl) in the full space. 
Secondly, there are forces that can be derived from a potential in the full space 
but not from a potential on the torus. We give examples of both cases in this 
section. 

4.1 Nonperiodic potential 

Consider again a particle diffusing on an infinite line, with potential U {x) = 
cos X. To this we add a constant force F: 

it = vt 
mvt = F + sina;t — ^mvt + \/2ni'^T 

The force F + sinx can be derived from the potential —Fx + cos a;, but this 
potential is not periodic. The relation ([s]) becomes here 

M{t) = ^ + - «o); - xo)) 'it^-^y^ - ^°)> 

Cl(t) C2(t) 

Note that I^F; (xt — xq)'^ — because is a constant. We have simulated 

this system for various values of the force F. In Fig. [3] the mobility, diffusion, 
the correlation of the displacement with the change in velocity (CI) and the 
correlation of the displacement with the integrated force (C2) are shown for 
F ~ 1.5 (a), F — 2.5 (b). In both cases the initial conditions are xo — vq — 
0. Note that all quantities converge to a constant value. This is because the 
diffusion and other correlations are defined as truncated correlation functions 
(see ([3])), so that the diverging parts of the correlations are subtracted. 

It is clear that the Sutherland-Einstein relation M = D/T is violated. Even 
though the correction term CI tends to zero for large times, the correction 
term C2 does not. However, the correction to the Sutherland-Einstein relation 
is not always equally big. E.g. for small forces the correction is small. For 
large driving as in Fig. |3] (b), the correction is also small. The force is then 
so big, that the potential U does not play an important role. Without the 
potential U the dynamics of the system is completely translation invariant, and 
with a simple change of coordinates one can prove that the Sutherland-Einstein 
relation holds. Therefore one expects that in the limit F — >■ 00 the mobility and 
diffusion are proportional, and have the same values as those for pure diffusion. 
This can indeed be seen in Fig. |4] (a) for the mobility and diffusion in a range 
of forces going from —5 to 5. For a force that is five times the amplitude of the 
potential, the mobility and diffusion are very close to 1, the expected values for 
a pure diffusion. Finally, as also visible in the figure, both M{t) and D{t) are 
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Figure 3: Simulations of a diffusion with a potential —Fx + cos a;. Here T = 
J — m = 1. In (a) = 1.5 and in (b) F = 2.5. (Online version in color.) 



symmetric under the force flip F — >■ —F. The resulting curves for mobility and 
diffusion are very similar to the overdamped case, as e.g. in Figure 3 of [22j . 
where the stationary probability density is exactly known. 




-4 -2 2 4 -8 -6 -4 -2 2 4 6 8 

F A 

Figure 4: The long time limits of the mobility and diffusion are shown. In (a): 
simulations of a diffusion with a potential —Fx + cos a;, where F ranges from 
—5 to 5. In (b) the potential is —A{x — cosx), with A ranging from —8 to 8. In 
both cases T = j — m = 1. (Online version in color.) 

Fig. |4] (b) shows the mobility and diffusion (long time limits) for a force 
A{1 + sin a:), with A between —8 and 8. Again, for large amplitude A the rela- 
tion approches that of pure diffusion just as for F — >■ oo in Fig. |4j The reason 
for that is unclear; the forcing is large but also the potential is large now. At 
any rate, the diffusion is clearly much more sensitive to the strength of the force 
than is the mobility. 
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Figure 5: The long time limits of the mobility and diffusion are shown for a 
potential given by —Fx + cosKx, with K ranging from to 10. Again, we took 
T — J — m — I. (Online version in color.) 



Finally, Fig. [s] shows the mobility and diffusion (long time limits) for a 
potential —Fx + cosKx, with K between and 10. For K between 1 and 10 
the difference between mobility and diffusion is almost constant. 

4.2 Rotational force 

Apart from forces coming from nonperiodic potentials, nonequilibrium can also 
be installed from forces with a nonzero curl on the full space. They induce 
vortices and let the particles undertake rotational motion, but with more drastic 



departures from the Sutherland-Einstein relation than in the case (17) of a 
magnetic field. 

As an example we simulated the two-dimensional dynamics 

mva:,t = Fx[xt,yt) ~ IVx.t + \/2m-fT^^^t (18) 
mvy^t = Fy{xt,yt) - JVy^t + \/^rnrff ^y^t 

The force F — Af, with A some constant and with, for < a:, y < 1, 
Ux,y) = ir-V2){^-y) 
fy{x,y) = (r-V2)(a;-i) 

with r — ^ {x — |)2 + (t/ — |)2. This is shown in the left part of Fig. 

This force field is repeated outside that unit square, however sometimes in 
the reverse direction. More precisely, 

fx{x,y) = a(r- \/2) (y- i) 

fy{x,y) = a{r~^/2){]^~x) 

a = (1 — 2S2^x mod 3) (1 ^ 2(5i J, mod 2) 

This gives the pattern shown in the right part of Fig. |6] The specific choices 
are probably not so important. Simulation results for this system are shown in 
Fig.[7l 
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Figure 6: Vector plot of the force F in (18 1. (Online version in color 
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Figure 7: Values T = 7 = m = 1. In (a) the constant A is 20, while in (b) 
A — 40. Initial conditions in both cases take position and velocity equal to zero. 
(Online version in color.) 



We only show the diagonal elements of the mobility, diffusion and the cor- 
rection (C); the off-diagonal elements turn out to be zero. In both cases the 
diffusion in the x-direction is bigger than in the y-direction. Furthermore, when 
comparing v4 = 20 with A = 40, we see that the diffusion very much depends on 
the strength of the force, while the mobility is in both cases approximately the 
same. To see this more clearly, Fig. [8] shows the large time limits of mobility, 
diffusion, and their difference, for A in the range [—60, -1-60]. 

Observe that the diffusion increases rapidly with the amplitude \A\ of the 
force, while the mobility remains almost constant (it even decreases a little). 
In particular, large forcing A f 00 does not at all lead to pure diffusion here, 
because the situation is very different from that in Fig. |4] We also see that all 
quantities are symmetric under A — )■ —A. Indeed, reversing the force is the same 
as translating the whole system, which has no effect on the long time behavior. 
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y = m = 1, and A ranges from —60 to +60. (Online 

5 The symmetrized mobility 

We now consider the symmetric part of the mobility matrix, i.e., we rewrite 
equation ([5| to: 

2 J. 4^^y^ 



(ft-fo),;^'j-) + ((r=i-fo)j;«' 



The last term can be rewritten using elementary algebra 
— ({rt- ro),; + —(irt- ro)^; * 

7717 ^ ' ^ 

= ( (?^t - H ; in- ro) 

\ L my J L 



in - ro)t; in - ro)] 



777^7^ 



7717 



The second term on the right-hand side is proportional to the diffusion function. 
Furthermore, the first of the terms on the right-hand side of the last equation 
can be simplified: for that, let us write the Langevin equation in its integral 
form: 

-F{rs,Vs) 



Vt-vo^ -lin - rj)) + / ds 
Jo 

Using definition Q of the functional 5*, wc substitute 



27T 



6 



-'i' + ft-fo 




ds 



The right-hand side involves Gaussian white noise and (^S.u.iCs,j^ — S{s — u)6i j. 
The result is the following expression for the symmetric part of the mobility: 



2 2T 27777 47277i2Tt\ " ^ 



(19) 
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Using exponential relaxation in large time limits (see below in Appendix IaI) one 



can see that in the long time limit this relation becomes ( 10 ). However, even for 
finite times when i = j, the third term on the right-hand side in ( |19[ ) is minus 
a variance, yielding the following bound for the mobility (diagonal elements): 

M,Jt) < + — (20) 

^ ' - 2T 2m7 ^ ' 

with equality for pure diffusion. 

6 Conclusions 

This work investigates the relation between diffusion and mobility for Langevin 
particles. The particles are independent and passive, undergoing white noise and 
friction from the fluid in thermal equilibrium. The total force on the particles 
is periodic but not compatible with a periodic potential, bringing the system 
out-of-equilibrium. We have studied the modified Sutherland-Einstein relation. 
The new mobility-diffusion relation remains explicit and we have visualized the 
typical dependencies on the nonequilibrium driving and other parameters. For 
the future we hope that the inverse analysis will also be possible, i.e., that 
our formula will enable to obtain useful information about unknown aspects of 
diffusive nonequilibria exactly by measuring the correction to the Sutherland- 
Einstein relation and from comparing it with ([t]) or with (10 1. In many cases 



however we expect that further extensions to nonMarkovian evolutions will be 
necessary to meet the physics of small systems immersed in viscoelastic media. 

A Smoothness and mixing 

We argue here why the diffusion matrix Dij (t) has a finite limit and why certain 
terms do not contribute in the long time limit; that is the vanishing of ^ 



and (15). We do not want to go into full technical details but it should be 
clear that our results depend on good exponential mixing properties of the 
dynamics with propagation of smooth densities. Of course, there is no stationary 
regime in unbounded diffusive systems, but this problem is irrelevant whenever 
expectations are considered of quantities that have the same periodicity as the 
dynamics (forces) of the system. (This gives no restriction on the velocity 
dependence of these quantities.) In that case we can restrict the dynamics 
of the system to one period, with periodic boundary conditions (reducing the 
infinite space to a torus). Taken on a torus, our system does have a stationary 
regime, i.e., the Fokker-Planck equation ^ with periodic boundary conditions 
can in principle be solved with left-hand side zero. The solution is assumed to 
be the smooth density p{r, v), giving the stationary distribution of positions (on 
the torus) and velocities. By assumption then, any function g with the same 
periodicity as the dynamics satisfies 

lim {g{ft,vt)) = / drdv p{f,v)g{f,v) = (g) 
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and there exist finite constants C and a > such that for large enough times 
t > 0, 



9 \ —at 

,5 ) e 



From the above bound, one can deduce bounds on truncated correlation func- 
tions. For example, writing g{t) = g{rti ^t), for t > s, 



where stands for the conditional expectation of g at time t, given the 

state at time s. We can then insert the exponential bounds: 



C 



1 + C'e- 



C 



p \ I p 

With this bound, one obtains that the diffusion matrix has a finite limit: 
lim \DiAt)\ < lim — 



(21) 
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Similarly, for ([9]), we find that 



< lim — 

t->00 2t Jq 



du 



goes to zero, by direct use of (21 1. 



The proof of equation (15) is similar to the one for ([9|, with one extra ele- 
ment, namely the fact that the dynamics in this case is an equilibrium dynamics, 
meaning that the associated stationary distribution p is the equilibrium distri- 
bution. Equilibrium is characterized by time-reversal invariance. This means in 
particular that the expectation value of a time-antisymmetric quantity is zero 
in equilibrium. In our case, this means that 



{{rt-ro)^; f dsVrVir)) =0 



The rest of the proof of (151 is quite straightforward. 
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B The mobility for one-dimensional diffusions 



We give a short proof of formula (16 1. Consider the stationary probability 
current for one-dimensional overdamped diffusion on the circle in the presence 



of a constant force /: 



= X[f - U']p~xTp 



where p is the stationary distribution. Because of stationarity, the current does 
not depend on the spatial coordinate x. We divide the last equation by p and 
integrate over the circle: 



dx 
p{x) 



Rxf 



Up to first order in /, the stationary current is thus given by 



Rxf 

dx 



I 



where Pq = e '^(^)/'^/ J e '^(^^/^ is the equilibrium distribution. This gives the 
mobility 



M = R — 



R'x 



because 



1 d f 
M = lim - — {xt - xq)-' 

t-^oo t of 



^ J g+U(x)/T J g~U(x)/T 

19/"* 
1™ 7 / (dxt) 



f=o 



C About the simulations 

All graphs in this paper are the result of simulations made with the programming 
language CH — h by directly applying the Langevin equation. This means (taking 
the simple example of one dimension) taking variables x and v which at each 
timestep of length dt change by the operations 

X — >■ X + vdt 

V ^ v + F{x, v)dt - ^vdt + ^j2m^TdtS^ 

where ^ is a random number drawn with a standard normal distribution (mean 
zero and variance 1). In this way a trajectory of consecutive positions and ve- 
locities is generated. In the same way trajectories are generated for a dynamics 
where a small constant force is added. The relevant quantities (diffusion, inte- 
grated forces, displacement of the particle) are then computed at each timestep 
by averaging over many simulated trajectories. The size of the timestep, of 
the force, and the number of simulated trajectories varied between the different 
examples in this text to find each time a good compromise between statistical 
accuracy and computation time. The length of the trajectories was taken such 
that the quantities showed a clear convergence to a constant value. As an ex- 
ample, for the case of the nonperiodic potential, we took dt — 0.002, / = 0.03, 
the number of timesteps 15000 and the number of trajectories 160000. As the 
goal of this paper is not providing accurate quantitative information, but rather 
a qualitative visualization of the relation between mobility and diffusion, we 
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did not compute the statistical and systematical errors. Rather, we check each 
time that the mobility equals the sum of the terms on the right-hand side of 
([5]). The mobility needs much more computation time (number of simulated 
trajectories) than the other quantities. The problem there is that one subtracts 
two quantities that are close to each other, and then divides by a small number. 
Furthermore, for the mobility one needs to simulate two systems: the perturbed 
and the unperturbed one. This also means that our results are numerically use- 
ful, as one does not need to compute the mobility separately. We used this 
fact in the simulations where the force was varied, as in Figures |4) [5) [8) The 
simulations there took a longer time, as each different value of the force required 
a different simulation. By leaving out the mobility, the number of trajectories 
could however be reduced (e.g. for the nonperiodic force the number was 8000). 
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